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ABSTRACT 

We discuss the spectrum arising from synchrotron emission by fast cooling (FC) 
electrons, when fresh electrons are continually accelerated by a strong blast wave, into 
a power law distribution of energies. The FC spectrum was so far described by four 
power law segments divided by three break frequencies v sa < v c < v m . This is valid 
for a homogeneous electron distribution. However, hot electrons are located right after 
the shock, while most electrons are farther down stream and have cooled. This spatial 
distribution changes the optically thick part of the spectrum, introducing a new break 
frequency, v ac < u sa , and a new spectral slope, F v oc i/ 11 / 8 for u ac < v < u sa . The 
familiar F v oc u 2 holds only for v < v ac . This ordering of the break frequencies is 
relevant for typical gamma-ray burst (GRB) afterglows in an ISM environment. Other 
possibilities arise for internal shocks or afterglows in dense circumstellar winds. We 
discuss possible implications of this spectrum for GRBs and their afterglows, in the 
context of the internal-external shock model. Observations of F u oc z/ 11 / 8 would enable 
us to probe scales much smaller than the typical size of the system, and constrain the 
amount of turbulent mixing behind the shock. 

Subject headings: radiation mechanisms: nonthermal — shock waves — gamma rays: 
bursts — turbulence 



1. Introduction. 



The spectrum of Gamma-Ray Bursts (GRBs) and their afterglows is well described by syn- 
chrotron and inverse Compton emission. It is better studied during the afterglow stage, where we 
have broad band observations. The observed behavior is in good agreement with the theory. Within 
the fireball model, both the GRB and its afterglow are due to the deceleration of a relativistic flow. 
The radiation is emitted by relativistic electrons within the shocked regions. According to the 
internal-external shock scenario, it has been shown (Fenimore et. al. 1996; Sari & Piran 1997) that 
in order to obtain a reasonable efficiency, the GRB itself must arise from internal shocks (ISs) within 
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the flow, while the afterglow is due to the external shock (ES) produced as the flow is decelerated 
upon collision with the ambient medium. In the simplest version of the fireball model, a spherical 
blast wave expands into a cold and homogeneous ambient medium (Waxman 1997; Meszaros, P. 
& Rees, M. 1997; Katz & Piran 1997; Sari, Piran & Narayan 1998, hereafter SPN). An important 
variation is a density profile p(r) oc r -2 , suitable for a massive star progenitor which is surrounded 
by its pre-explosion wind. 

In this letter we consider fast cooling (FC), where the electrons cool due to radiation losses on 
a time scale much shorter than the dynamical time of the system 3 , tdyn- Both the highly variable 
temporal structure of most bursts and the requirement of a reasonable radiative efficiency, suggest 
FC during the GRB itself (Sari, Narayan & Piran 1996). During the afterglow, FC lasts ~ 1 hour 
after the burst for an ISM surrounding (SPN; Granot, Piran Sz Sari 1999a) and ~ 1 day in a dense 
circumstellar wind environment (Chevalier & Li 1999). We assume that the electrons (initially) and 
the magnetic field (always) hold fractions e e and e# of the internal energy, respectively. We consider 
synchrotron emission of relativistic electrons which are accelerated by a strong blast wave into a 
power law energy distribution: N(^y) oc 7~ p for 7 > j m = (p — 2)e e e'/(p — l)n'm e c 2 = e e e' / 3n' m e c 2 , 
where n' and e' are the number density and internal energy density in the local frame and we have 
used the standard value, p = 2.5. 

After being accelerated by the shock, the electrons cool due to synchrotron radiation losses. 
An electron with a critical Lorentz factor, j c , cools on the dynamical time, td yn 

7 C = &TTm e c/a T TB' 2 tdyn = 3m e c/A(j T Y£Be'tdyn , (1) 

where V is the bulk Lorentz factor, <jt is the Thomson cross section and B' is the magnetic field. 
The Lorentz factors, 7 m and 7 C , correspond to the frequencies: v m and v c , respectively, using 
v S yn{l) = 3g e 5'7 2 r/16m e c. FC implies that 7 C <C 7 m and therefore v c <C v m . 

The FC spectrum had so far been investigated only for v sa < v c < v m , using a homogeneous 
distribution of electrons (SPN, Sari Sz Piran 1999), where v sa is the self absorption frequency. 
Under these assumptions the spectrum consists of four power law segments: F v oc u 2 , z^ 1 / 3 , v^ 1 / 2 
and z/~ p / 2 , from low to high frequencies. The spectral slope above v m is related to the electron 
injection distribution: the number of electrons with Lorentz factors ~ 7 is oc 7 1_p and their 
energy oc 7 2 ~ p . As these electrons cool, they deposit most of their energy into a frequency range 
~ v S yn{l) oc 7 2 and therefore F v oc j~ p oc v~ p l 2 . At v c < v < v m all the electrons in the system 
contribute, as they all cool on the dynamical time, td yn . Since the energy of an electron oc 7, 
and its typical frequency oc 7 2 the flux per unit frequency is oc 7 _1 oc v~ x l 2 . The synchrotron 
low frequency tail of the cooled electrons (oc zW 3 ) appears at v sa < v < u c . Below v sa , the the 
system is optically thick to self absorption and we see the Rayleigh-Jeans portion of the black body 
spectrum: 

Fv OC V 2 lty V {v) , (2) 



s t dyn is the time for considerable expansion. Adiabatic cooling is therefore negligible compared to radiative cooling. 
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where 7* 2/p (^) is the typical Lorentz factor of the electrons emitting at the observed frequency v. 
Assuming Jtyp(v) = 7 C oc u°, one obtains F u ecu 2 . 

We derive the FC spectrum of an inhomogeneous electron temperature distribution in §2 . We 
find a new self absorption regime where F u oc u 11 ^ 8 . In §3 we calculate the break frequencies and 
flux densities for ESs (afterglows) with a spherical adiabatic evolution, both for a homogeneous 
external medium and for a stellar wind environment. ISs are treated in §4. In section §5, we 
show that the early radio afterglow observations may be affected by the new spectra. We find that 
synchrotron self absorption is unlikely to produce the steep slopes observed in some bursts in the 
1 — 10 keV range. We also discuss the possibility of using the new spectra to probe very small 
scales behind the shock. 



2. Fast cooling spectrum 

The shape of the FC spectrum is determined by the relative ordering of v sa with respect to 
v c < v m . There are three possible cases. We begin with u sa < v c < v m , case 1 hereafter. This is 
the "canonical" situation, which arises for a reasonable choice of parameters for afterglows in an 
ISM environment. The optically thin part of the spectrum (y > u sa ) of an inhomogeneous electron 
distribution is similar to the homogeneous one. All the photons emitted in this regime escape the 
system, rendering the location of the emitting electrons unimportant. In the optically thick regime 
{y < i^sa), most of the escaping photons are emitted at an optical depth r u ~ 1, and ^typ{v) must 
be evaluated at the place where r v = 1. 

In an ongoing shock there is a continuous supply of newly accelerated electrons. These electrons 
are injected right behind the shock with Lorentz factors 7 > j m , and then begin to cool due to 
radiation losses. In the relativistic shock frame, the shocked fluid moves backwards at a speed of 
c/3: I' = ct'/3, where I' is the distance of a fluid element behind the shock and t' is the time since 
it passed the shock. Just behind the shock there is a thin layer where the electrons have not had 
sufficient time to cool significantly. Behind this thin layer there is a much wider layer of cooled 
electrons. All these electrons have approximately the same Lorentz factor: 7(t') = &irm e c/ aTB' 2 t' , 
or equivalently, 7(T) = 2-7rm e c 2 /<jtB ,2 1' . Electrons that were injected early on and have cooled 
down to 7 C are located at the back of the shell, at a distance of A' = ct' dyn /3 = cTtdy n /3 = TA 
behind the shock (i.e. 7 C = 7(A') = 7(i^„))- We define the boundary between the two layers, 1' , 
as the place where an electron with an initial Lorentz factor 7 m cools down to 7 m /2 

l' = 2Trm e c 2 /ajnB ,2r y m = 3m 2 c 4 n' / AaT^e^B^' 2 ■ (3) 

The un-cooled layer is indeed very thin, as l' /A' = jc/lm *C 1. We define by t u (1[) = 1. The 
optically thin emission from /' < l[ equals the optically thick emission: 

nAT 2 P v , max [v/vsyMl^))} 1/3 /4tt = {2u 2 /c 2 ) r 7 (/i)m e c 2 , (4) 
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where P Vjma x ~ Psyn{l) /^syn{l) oc 7° and P syn {l) = r 2 07C7 2 -B /2 /67r are the peak spectral power 
and total synchrotron power of an electron, respectively. Since v S yn(l) oc 7 2 > and within the cooled 
layer, 7 = oc l/l', eq. (4) implies jt yp = oc v~^^ ■ We now use eq. (2) and obtain that 
F v oc i/ 11 / 8 . This new spectral regime, is a black body spectrum, modified by the fact that the 
effective temperature i^ftyp) varies with frequency. 

At sufficiently low frequencies, < 1' , implying ^ t yp = 7m oc u° and F u oc v . The transition 
from absorption by the cooled electrons (F v oc z/ 11 / 8 ) to absorption by un-cooled electrons (F u oc v 2 ) 
is at u ac , which satisfies V^Uac) = ^o- The resulting spectrum is shown in the upper frame of Fig. 1. 
v ac may be obtained from eq. (4) by substituting 7(7' 1 ) = 7m and l[ = 1' from eq. (3). Substituting 
7(^1) = 7c from eq. (1) and l[ = A' = TA into eq. (4), gives us v^} ■ The superscript, W, labels 
the specific case under consideration. We obtain 

v ac - [pm e c 7 n /7T e B e e q e e ) , 

ug) = (8/3tt) [44e|r 13 AV 3 e' 6 /9m^c 3 ?e 4 ] 1/5 . (5) 

The ratio v)t} jv ac = (7 m /7c) 8 ^ 5 = (^m/^c) 4 ^ 5 depends on the cooling rate. The maximal flux 
density occurs at v c = v S yn{lc) and is given by F^max = N e P u>max /A7rD 2 (SPN), where N e is the 
number of emitting electrons and D is the distance to the observer, while v m = v syn (^ m ) 

u c = 27^m e cq e / UV2a 2 T Te% 2 t 2 dyn e^ 2 , 

Um = ^q e Te]i 2 e 2 e e' 5 / 2 / 12V2m 3 e c 5 n' 2 , (6) 
Fornax = ^V2a T m e c 2 N el e]i 2 e n / 2 / ^' 2 q e D 2 . 

For v c < u sa < v m (case 2) the cooling frequency, u c , becomes unimportant, as it lies in the 
optically thick regime. Now there are only three transition frequencies. u ac and v m are similar to 
case 1. The peak flux, F„% ax is reached at via '■ 

= (^ } ) 5/ V /9 , 
F^ = F« (v /u^Y /18 (7) 

-*■ u,max 1 u,max I "c/ "so I \ ' ) 

If v c < v m < v sa (case 3) then l[ (v) <C 1' for v < . Now, both u ac and v c are irrelevant, as the 
inner parts, where these frequencies are important, are not visible. We can use the initial electron 
distribution to estimate Jtyp- Jtyp = 7m °c v° at v < v m , implying F u oc v 2 . At u rn < v < v sa the 
emission is dominated by electrons with v syn (7) ~ 17 implying j typ oc u 1 / 2 and F u oc v^l 2 . F^)n ax 
is reached at vf}'- 



„(3) 

sa 



\"sa ) c m 



i?(3) = F (l) 



(^))-/3,(3-P)/ 6 ^-l)/2l 5/(P+5) . ( 8) 
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3. Application to External Shocks and the Afterglow 

Consider now the FC spectrum of an ES which is formed when a relativistic flow decelerates 
as it sweeps the ambient medium. This is the leading scenario GRB afterglow. We consider an 
adiabatic spherical outflow running into a cold ambient medium with a density profile p(r) oc r~ a , 
for either a = (homogeneous ISM) or a = 2 (stellar wind environment). FC lasts for the first hour 
or so in a typical ISM surrounding, and for about a day in a stellar wind of a massive progenitor. 
The proper number density and internal energy density behind the shock are given by the shock 
jump conditions: n' = 4Tn and e' = 4T 2 nm p c 2 , where n is the proper number density before the 
shock and m p is the mass of a proton. We also use A = R/12T 2 , R = 4F 2 ct and td yn = t, where t 
is the observed time. 

For a homogeneous environment iV e = 4irnR 3 /3. Using T oc i~ 3 / 8 (e.g. SPN) we obtain 

- 17v1(l9R 7 -2/5 -8/5^-1/10 3/10 3/10 



„£> = 1.8 x 10- Hz nl^V 72 , 



2.9 x 10 15 Hz e^E^n^H, 



2 

_ = = x ln 16 tt 7 1/2 2 #1/2,-3/2 



F^l ax = 30 mJy D^ef 01 E 52 n\ /2 , (9) 

where n\ = n/1 cm" 3 , E 52 = E/10 52 ergs, D 28 = D/W 28 cm, e.B,o.i = es/0.1, e e ,o.i = e e /0.1 and 
t 2 = t/WO sec. For typical parameters, only the case 1 spectrum is expected. After ~ 1 hour, slow 
cooling (SC) sets in, and the spectrum is given in GPS. 

For a circumstellar wind environment, n = r~ 2 A/m p and N e = 4:irAR/m p . Using T oc i -1 / 4 
and A = 5 x lO 11 ^ gr cm -1 as in Chevalier & Li (1999) we obtain 

V - 3 6 X 10 10 Hz 4 3 /5- 2 /5 -8/5 „-2/5 
Vac — o.O X 1U JTLZ i-B,Q.\ e efS.\ ^52 > 

- 8 x 10 11 Hz4 11/5 e 6/5 E~ 2/5 t~ 8/5 

v sa — °- u x iu nzA * e S,0.1- C/ 52 E /ir ' 

!/ c = 2.5xl0 12 Hz4- 2 e -j4 /2 if , 
*/ m = 1.2 x 10 14 Hz e^'o.i^C 72 , 

- n T,r D-2 4 C V2 pl/2 -1/2 
r i/,max — u - oy J J ^28 ■ f± * t B,0.1 I2j 52 L hr ' 

i£> = 1.3 x 10 12 Hz^/ 3 t- 2/S > 

= 0.54 Jy D^-^/^tV" , (10) 

where ^ r = t/1 hour. For typical parameters, the spectrum is of case 2 for ~ 1 — 2 hours after the 
burst. Then it turns to case 1 until ~ 1 day, when there is a transition to SC. The SC spectrum is 
given in Chevalier & Li (1999). 
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4. Application to Internal Shocks and the GRB 

ISs are believed to produce the GRBs themselves. The temporal variability of the bursts is 
attributed to emission from many different collisions between shells within the flow. The number 
of peaks in a burst, N, roughly corresponds to the number of such shells. Different shells typically 
collide before their initial width, Aj, has expanded significantly. Assuming that the typical initial 
separation between shells is ~ Aj, Aj = cTgo/2N in average, where T90 is the duration of the burst. 
The average energy of a shell is E s h ~ E/N, where E is the total energy of the relativistic flow. The 
emission in the optically thick regime comes from the shocked fluid of the outer and slower shells. 
We denote the initial Lorentz factor of this shell by Tj, and its Lorentz factor after the passage of 
the shock by V. The average thermal Lorentz factor of the protons in this region equals the relative 
bulk Lorentz factor of the shocked and un-shocked portions of the outer shell, r r = r/2Tj, which 
is typically of order unity. Therefore e' = r r n'm p c 2 . One can estimate the number density of the 
pre-shocked fluid, n^, by the number of electrons in the shell, N e = Egh/TimpC 2 , divided by its 
volume: = N e /4irR 2 AiTi. The number density of the shocked fluid, which is the one relevant 
for our calculations, is n' = 4r r n'j. The width of the front shell in the observer frame decreases 
after it is shocked: A = Aj/8r 2 . In this section we use R ~ 2r?A« ~ 4r 2 A, which is the typical 
radius for collision between shells, and td yn ~ 3A/c. Thus, we obtain 

„ ac = 8.7 x 10- Hz r^^^^^r^ivfT; 3 / 5 , 
i&> = 7.7 x io» Hz rf/ 5 41^3 41/5 iv|rr 19/5 , 

v c = 3.6 x 10 12 Hz V;y B Z l\E-^ 2 TlN^Tl /2 , 
v m = 5.6 x 10 18 Hz Tf^l^E^TfN^T-^ 2 , 
FlZLx = 509 Jy ry^EH 2 rj 3 £>^ 2 T^' 2 , 

v?) = 4.3 x 10 16 Hz r^El^N^Tr 1 , (11) 
where N 2 = AT/100, Ti = T 90 /10 sec, r rj3 = T r /3 and T 3 = T/10 3 . 



5. Discussion 

We have calculated the synchrotron spectrum of fast cooling (FC) electrons. We find three 
possible spectra, depending on the relative ordering v sa with respect to v c < v m . Two of these 
spectra contain a new self absorption regime where F v oc z/ 11 / 8 . 

During the initial fast cooling stage of the afterglow, the system is typically optically thick in 
the radio and optically thin in the optical and X-ray, for both ISM and stellar wind environments. 
We therefore expect the new feature, F v oc z^ 11 / 8 , to be observable only in the radio band, during 
the afterglow. For both environments, v ac and v sa move closer together (see Fig. 1) until they 
merge at the transition to slow cooling (SC). Afterwards, there is only one self absorption break, 

at fsa ^ ' 
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For an ISM surrounding, u sa oc t . From current late time radio observations, we know that 

( SO 

typically, u sa ~ a few GHz (Taylor et. al. 1998; Wijers & Galama 1999; Granot, Piran Sz Sari 
1999b). Therefore, the whole VLA band, 1.4 — 15 GHz, should initially be in the range where 
F v oc z/ 11 / 8 . Sufficiently early radio observations, which could confirm this new spectral slope, may 
become available in the upcoming HETE era. 

In a considerable fraction of bursts, there is evidence for a spectral slope > 1/3 (photon number 
slope > —2/3), in the 1 — 10 keV range (Preece et. al. 1998; Crider et. al. 1998; Strohmayer et. al. 
1998). Such spectral slopes are not possible for optically thin synchrotron emission (Katz 1994). 
They could be explained by self absorbed synchrotron emission if v sa reaches the X-ray band. The 
best prospects for this to occur are with the spectrum of the second type. However, we have to 
check whether the physical parameters for which v sa is so high are reasonable. Several constraints 
must be satisfied: (i) ISs must occur at smaller radii than ESs, (ii) efficient emission requires FC, 
and (iii) the system must be optically thin to Thomson scattering and pair production. The most 
severe constraint in the way of getting the v sa into the X-ray band arises from (iii) . It is possible 
only for rather extreme parameters: T ~ 10 4 and T r ~ 50. With such parameters, vF v would peak 
at hv, m ~ 1 — 100 GeV, unless e e < 10~ 2 , which would result in a very low radiative efficiency. 
Overall it seems unlikely that self absorbed synchrotron emission produces the observed spectral 
slopes > 1/3. 

So far we have neglected Inverse Compton (IC) scattering. For e e < e# the effects of IC 
are small, as the total power, Pic, emitted via IC is smaller than via synchrotron: Pic < P S yn- 
For e e > eb IC becomes important as Pic/Psyn = \AeA-B (Sari, Narayan & Piran 1996). This 
additional cooling causes 7 C and v c to decrease by factors of \Je e j^B and e e /es, respectively. This 
increases the duration of the FC stage in ESs by a factors of e e /ee and y 'e e /eB for ISM and stellar 
wind environments, respectively. 

Our results depend on the assumption of an orderly layered structure behind the shock: a 
thin un-cooled layer of width 1' , followed by a much wider layer of cooled electrons, of width 
A' = ZgyV m /f c . Clearly, significant mixing would homogenize the region and would lead to the 
"homogeneous" spectrum given in SPN, without the zv 11 / 8 region discussed here. A typical electron 
is not expected to travel much farther than its gyration radius, r(^). We obtain r(7 m )/Zg = 
5 x lO" 8 ^ 9/8 , 10~ 9 t^ r 5/4 and 5 x 10~ 8 , using the scalings of eqs. (9) and (10) for ESs and eq. (11) 
for ISs, respectively. Thus, this effect could not cause significant mixing. Another mechanism that 
might cause mixing is turbulence. The observation of a spectrum with F v oc v 11 / 8 would indicate 
the existence of a layered structure, and constrain the effective mixing length to < 1' = IS! ^Jv c jv m . 
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Fig. 1. — Fast cooling (FC) synchrotron spectra from a shock injected power law electron distri- 
bution. The shape of the spectrum is determined by the ordering of the self absorption frequency, 
u sa , with respect to u c < v m . There are three possible shapes for the spectrum, corresponding to 
v sa < v c , u c < v sa < u m and v sa > v m , from top to bottom. Scalings for an afterglow, both in an 
ISM and wind environments, are given only in the first two frames, since for typical parameters 
v sa < v m during the FC phase of the afterglow. 



